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ABSTRACT 

The redshift distribution of near-IR selected galaxies is often used to attempt to discriminate between 
the classical view of galaxy formation, in which present-day luminous galaxies were assembled at early 
times and evolve due to the passive aging of their stellar populations, and that of hierarchical structure 
formation, in which galaxies were assembled more recently via the merging of smaller objects. We carry 
out such a test here, by computing the distribution of photometric redshifts of A"ab < 22 galaxies in the 
Great Observatories Origins Deep Field Survey (GOODS) Southern Field, and comparing the results 
with predictions from a semi-analytic model based on hierarchical structure formation, and a classical 
'passive evolution' model. We find that the redshift distributions at z < 1.5 of both the hierarchical 
and passive models are very similar to the observed one. At z > 1.5, the hierarchical model shows a 
deficit of galaxies, while the passive model predicts an excess. We investigate the nature of the observed 
galaxies in the redshift range where the models diverge, and find that the majority have highly disturbed 
morphologies, suggesting that they may be merger-induced starbursts. While the hierarchical model used 
here does not produce these objects in great enough numbers, the appearance of this population is clearly 
in better qualitative agreement with the hierarchical picture than with the classical passive evolution 
scenario. We conclude that the observations support the general framework of hierarchical formation, 
but suggest the need for new or modified physics in the models. 

Subject headings: galaxies: formation — galaxies: evolution 



O ,■ 1. INTRODUCTION 

q While the stars that produce the bulk of the optical light 
?-h ' in nearby luminous galaxies are known observationally to 
c/3 , be quite old, it is possible that the time at which the mass 
C3 ■ was assembled into a single object is different from the 
£^ | formation time of the stars. In fact, in the hierarchical 
• i-H . paradigm of galaxy formation, one expects mass assembly 
' to be a gradual process, in contrast to the classical mono- 
H | lithic dissipative collapse picture (Eggen, Lynden-Bell, & 
. - - i Sandage 1962). A generic prediction of hierarchical (Cold 
Dark Matter; CDM) theories is that galaxies should be less 
massive in the past. Direct searches for galaxies at high 
redshift should therefore provide a crucial test of this class 
of theories (Kauffmann & Chariot 1998). Determining the 
epoch of formation and assembly of present-day luminous 
galaxies is central to our understanding of galaxy forma- 
tion and cosmology, and is a primary goal of the Great 
Observatories Origins Deep Survey (GOODS), as well as 



of many other deep surveys. 

In practice, however, there are many complications in- 
volved in carrying out this test. Of course, we cannot fol- 
low an individual galaxy back in time, but can only study 
populations observed at different redshifts. Cosmological 
k-corrections (redshifting of light to longer wavelengths) 
and stellar evolution both make it difficult to relate high 
redshift populations to local ones. One can attempt to ac- 
count for these effects by artificially redshifting template 
spectra for representative galaxy types and/or by using 
stellar population models with assumed star formation his- 
tories. Starting from an empirical low-redshift luminosity 
function, one can then predict how the local galaxy pop- 
ulation would appear at high redshift, with the effects of 
k-corrections only ('no evolution') or with the additional 
effects of stellar evolution ('passive evolution') included. 
However, these corrections can be sensitive to the input 
assumptions. 

There are several advantages to carrying out this test 
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using a sample selected in the near-IR rather than the 
optical. The corrections described above are consider- 
ably smaller, and the near-IR light more closely traces 
the stellar mass. The observed K-band (2.2 /zm) probes 
the SED longwards of the rest-frame I-band (<~ 8600 A) 
out to z ~ 1.5. However, it also necessary to probe a 
large enough volume to contain a statistically significant 
sample of rare luminous objects. The availability of accu- 
rate photometric redshift estimates, requiring multi-band 
U through K photometry, enables such tests to be ex- 
tended into the z > 1.2 'spectroscopic desert', where spec- 
troscopic redshifts are difficult to procure. Only recently 
have sufficiently deep and wide near-IR selected surveys 
with multi-band photometry begun to become available. 

The original study by Kauffmann & Chariot (1998, 
KC98) showed that the cumulative redshift distribution of 
K-selected galaxies differed greatly in a hierarchical model 
compared with a passive (or Pure Luminosity Evolution; 
PLE) model: they found that the fraction of galaxies at 
z > 1 with 18 < i^vEGA < 19 was an order of magnitude 
higher in the PLE model. Based on the small observa- 
tional samples available at the time, KC98 concluded that 
the hierarchical model provided a better match to the ob- 
served redshift distribution out toz~ 1 . Recently, several 
studies have carried out this test using updated hierar- 
chical models (including, among other things, the transi- 
tion to a low-f2 m , cosmological constant-dominated cos- 
mology) and larger, deeper observational samples (Firth 
et al. 2002; Cimatti et al. 2002b; Kashikawa et al. 2003). 
Both Cimatti et al. (2002b) and Kashikawa et al. (2003) 
concluded that the PLE models produced better agree- 
ment with the observed redshift distribution at z > 1.5 
than the hierarchical models, although the disagreement 
was relatively subtle compared with the expectations set 
out in KC98. 

In this paper, we use a K s -band selected sample 
{K AB < 22) from the GOODS Southern field to repeat the 
KC98-type redshift distribution test, using accurate, well- 
calibrated photometric redshifts (Mobasher et al. 2003). 
The GOODS field probes a considerably larger area and 
volume than previous studies at a similar depth. We also 
have the advantage of the exquisite ACS imaging, which 
allows us to investigate the morphologies of high redshift 
galaxies, gaining further insights into the nature of these 
objects. We confront these observations with predictions 
from semi-analytic hierarchical galaxy formation models, 
and with PLE models, both normalized to the z = K- 
band luminosity function recently determined from the 
2MASS survey (Kochanek et al. 2001; Cole et al. 2001). 

2. THE DATA 

The GOODS data are described in Giavalisco et al. 
(2003). Our study is based on the GOODS Southern field 
(Chandra Deep Field South; CDFS), which has an area 
of 160 arcmin 2 . In addition to the 4-band (BViz) ACS 
imaging, we make use of an extensive set of complemen- 
tary ground-based observations from the VLT, NTT, and 
ESO 2.2m telescopes, including optical WFI (U'UBVRI) 
and FORS (RI), and infrared SOFI (JHK S ) photometry, 
which covers the entire ACS GOODS CDFS field. This 
work is based on a PSF-matched SOFI-K s selected catalog. 
Clearly-unresolved sources (stars) based on the ACS zsso 



data down to zab < 26.2 have been removed. The 50% 
completeness limit is K s = 22.8, and the sample should 
be close to 100% complete at K s < 22 (Moy et al. 2003; 
Giavalisco ct al. 2003), the limit we adopt for our analy- 
sis. Photometric redshifts were estimated as described in 
Mobasher et al. (2003), using all the available bands (IP 
through K s ), and are well-calibrated to i^vEGA < 20 using 
spectroscopic redshifts from the K20 survey (Cimatti et al. 
2002a) and additional spectra obtained from FORS2 on 
the VLT as part of the GOODS program. Typical redshift 
errors for the K s < 22 sample are ~ Az/(1 + z spcc ) = 0.1 
(Mobasher et al. 2003). 

In the remainder of the paper, we refer to the K s band 
as K for brevity, and give all magnitudes in the AB sys- 
tem unless otherwise specified. Note that for our fil- 
ter bands, K A b = Kvega + 1-85, and (R — K)ab = 
(R-K) VEGA -IM. 



3. MODELS 

Where relevant, we assume the following values for 
the cosmological parameters: matter density fl m = 0.3, 
baryon density Qb — 0.044, dark energy Cl\ = 0.70, Hub- 
ble parameter Ho = 70kms~ 1 Mpc -1 , fluctuation ampli- 
tude as = 0.9, and a scale-free primordial power spec- 
trum n s = 1. These values are consistent with the recent 
WMAP data (Spergel et al. 2003). 

The basic ingredients of the semi-analytic hierarchical 
models used here are described in Somerville & Primack 
(1999) and Somerville et al. (2001). The models are based 
on hierarchical merger trees within a ACDM model, and 
include modeling of gas cooling, star formation, supernova 
feedback, chemical enrichment, stellar population synthe- 
sis, and dust. We use the multi-metallicity stellar SED 
models of Devriendt et al. (1999), assuming a Kennicutt 
IMF. Here we have considered a model based on the 'col- 
lisional starburst' recipe described in Somerville et al. 
(2001), which was found to produce the best agreement 
with high redshift [z ~ 3) galaxy observations. Several 
parameters and model ingredients have been adjusted to 
give better agreement with the low redshift optical and K- 
band luminosity functions recently determined by SDSS 
and 2MASS (e.g. Blanton et al. 2001; Cole et al. 2001), 
and with low redshift galaxy colors (details will be given in 
Somerville et al. 2003, in prep). We produced a mock cata- 
log with the same angular extent and depth as the GOODS 
ACS and ground-based data, which was run 'blind' before 
the data were analyzed. 

The passive evolution models are computed as de- 
scribed in Gardner (1998), and are normalized to the type- 
dependent z = K-band luminosity functions derived 
from the 2MASS survey by Kochanek et al. (2001). These 
models contain six different types of galaxies with simple 
parameterized star formation histories: E, SO, Sb, Sc, Irr, 
and starburst. All galaxies except the starburst type be- 
gin forming stars at z = 15. The E, SO, Sb and Sc types 
have exponentially declining star formation rates with e- 
folding timescales of 1 Gyr for the E/S0, 4 Gyr for the 
Sb and 7 Gyr for the Sc types. The Irr types have a con- 
stant star formation rate. The starburst population has a 
constant star formation rate and constant age (1 Gyr) at 
every redshift. 
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Fig. 1.— [left] Cumulative redshift distribution of K < 22 galaxies. Bold solid line: GOODS CDFS; Light solid line: K20 sample; 
dashed: semi-analytic model; dot-dashed: passive model, [right] Differential redshift distribution (K < 22), with Poisson errors 
shown on the GOODS data [bold histogram]. Line types are as in the left panel. The hierarchical model is highlighted by diamond 
symbols. The cross symbols show the hierarchical model with a crude correction of 0.2 magnitudes applied to account for the 
isophotal magnitudes used in the GOODS observations (see text). 



4. RESULTS 

The cumulative and differential redshift distributions for 
the K < 22 selected samples are shown in Fig. 1, for 
the GOODS data, the K20 survey, and the semi-analytic 
and passive models. The redshift distribution of GOODS 
agrees very well with that obtained from the K20 spec- 
troscopic survey, which was carried out independently in a 
smaller area of the same field. The agreement between the 
semi-analytic models and the observations is quite good 
— well within the fluctuations expected from large scale 
structure — up to about z ~ 1.5. This good low red- 
shift agreement is in contrast with previous comparisons 
with semi-analytic models from several groups (Firth et al. 
2002; Cimatti et al. 2002b). In previous models, the lumi- 
nosity function was too steep on the faint end, leading to 
an excess of intrinsically faint galaxies at low redshift. Im- 
proved modeling of sub-halo merging, and including ejec- 
tion of gas by superwinds and suppression of gas infall 
in small halos after reionization produce better agreement 
with the observed K-band luminosity function at z = 
(Somerville 2002; Benson et al. 2002), and also to better 
agreement with the low-redshift N(z) shown here. Also of 
note is the similarity of the predicted N(z) at z < 1 for 
the hierarchical and the PLE models — in contrast with 
the results of KC98. 

At z > 1.5, the number of K < 22 galaxies in the semi- 
analytic models is significantly and systematically smaller 
than the observed value, while the PLE model systemat- 
ically overpredicts the number of objects in this range by 
a similar factor. The isophotal magnitudes used in our 
GOODS catalog are probably fainter than the true total 
magnitudes by about 0.2-0.3 magnitudes (Cimatti et al. 



2002a). Correcting the semi-analytic models for this effect 
further exacerbates the discrepancy, as shown in Fig. 1. 
We note that the redshift range in which the models suffer 
from the most significant discrepancy (z > 1.5) is precisely 
where the photometric redshifts are the least secure, as 
very few spectroscopic redshifts are available to test them 
— though progress is being made in this area (Daddi et al. 
2003). In addition, at these redshifts, the observed K-band 
has shifted into the rest optical, and is therefore more sen- 
sitive to recent star formation activity. 

It is interesting to investigate the nature of the galax- 
ies that appear in the high redshift tail of the observed 
distribution, but are under-represented in the hierarchi- 
cal models. Fig. 2 shows the (observed frame) R — K 
colors as a function of redshift for the observed and semi- 
analytic model galaxies. The colors of solar metallicity, 
single-age populations are also shown. The semi-analytic 
model reproduces the locus of observed colors fairly well 
up to z ~ 0.8. At z ~ 0.8-1.2, the semi-analytic model 
produces enough K < 22 galaxies overall, but does not 
produce enough galaxies with very red colors. This prob- 
lem has been noted before (Daddi et al. 2000; Firth et al. 
2002; Cimatti et al. 2002b). The number densities and 
morphologies of these R — K > 3.35 Extremely Red Ob- 
jects (EROs) in GOODS are discussed in more detail in 
Moustakas et al. (2003). 

At higher redshift (z > 1.5), the observed distribution of 
R — K colors is bimodal, with the 'dip' around R — K ~ 3 
(see Fig. 2). Focusing on galaxies in the redshift interval 
1.7 < z < 2.5 and with K < 22, we find that about 60% 
(44 out of 75) of the GOODS galaxies have colors bluer 
than R — K = 3, while in the semi-analytic mock catalog, 
52% (11 out of 21) of the objects have R - K < 3. Con- 
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sidering the small number statistics, this implies that the 
semi-analytic model produces approximately the correct 
relative fraction of red and blue galaxies. We have visu- 
ally inspected all of the objects in this sub-sample, and 
find that the great majority of both red and blue galaxies 
have highly irregular morphologies, many with multiple 
components and the appearance of ongoing mergers. 

From comparison with the single-burst model tracks, we 
can deduce that the high redshift, blue galaxies (R — K < 
3) must be dominated by extremely young (< 500 Myr), 
nearly unreddened stars. Intriguingly, Daddi et al. (2003) 
have recently obtained spectra for a sample of galax- 
ies with JsTvega < 20 and z p hot > 1-7 in the GOODS 
CDFS/K20 field, and have successfully obtained redshifts 
for 9 such objects, confirming that they lie in the range 
1-7 < z < 2.3. On the basis of these spectra and the ACS 
images, Daddi et al. (2003) argue that these objects are 
strongly clustered, massive, merger-driven starbursts. We 
show the location of these objects on our color-redshift di- 
agram, and see that they lie precisely in the regime of the 
'missing blue galaxies'. 

5. DISCUSSION 

In this paper, we address a central question in galaxy 
formation theory: were most of the luminous galaxies that 
we see today already in place at high redshift, or were they 
assembled gradually over time? To answer this question, 
we confronted observations from the GOODS CDFS field 
with predictions from two models representing what are 
traditionally considered opposing points of view: a semi- 
analytic hierarchical model based on CDM theory, and 
a 'passive evolution' model, in which galaxy properties 
evolve only due to the aging of their stellar populations. 
Our main conclusions are as follows: 

1. Up to z ~ 1, redshift distributions of K < 22 galax- 
ies in the hierarchical model, the passive model, 
and the data are all consistent with one another. 
However, the hierarchical model underproduces the 
number of Extremely Red Objects (EROs) at z <~ 1. 

2. The hierarchical model underproduces near-IR se- 
lected objects (K < 22) by about a factor of three at 
z > 1.7 and by an order of magnitude at z > 2. The 
PLE model overproduces these galaxies by about a 
factor of two at z ~ 2. 

3. At z > 1.5, the objects underproduced in the hier- 
archical model are nearly equally divided between 
red (R — K > 3) and blue galaxies. Based on ACS 
imaging, many of these objects appear to be highly 
morphologically disturbed, and a large fraction may 
be merger-driven starbursts. 

Not surprisingly, the predicted colors of model galaxies 
in the hierarchical models are quite sensitive to the details 
of the star formation recipes, as well as the stellar IMF 
and dust modeling. For example, if we assume that star- 
bursts occur in major mergers only, we can produce more 



extremely red galaxies at z ~ 1, but we then produce even 
fewer luminous blue galaxies at z ~ 1.5-2. Alternatively, 
if we brighten all model galaxies by 0.5 magnitudes (40%), 
we find that the semi-analytic model then produces suffi- 
cient numbers of K < 22 objects at z > 1.5, but this natu- 
rally causes an excess at lower redshift. Cosmic variance is 
also expected to be significant for these luminous, rare ob- 
jects — assuming that these objects are strongly clustered, 
like EROs at z ~ 1 (e.g. Daddi et al. 2001), we estimate 
an uncertainty due to cosmic variance of about 60% in the 
number density of objects at 1.5 < z < 2 (see Somerville 
et al. 2003). This implies that the semi-analytic model is 
discrepant at less than 2a. Results from additional fields 
will determine whether there is an overdensity of objects 
at z > 1.5-2 in the CDFS. 

A significant conclusion from this work is that the test 
proposed by KC98, when carried out with recent models, 
is not as strong a discriminator between the traditionally 
opposing points of view of hierarchical vs. PLE models as 
was found in that work. Adoption of the flat, \ow-Q m cos- 
mology now favored by observation, and the refinement 
of the star formation and feedback recipes has resulted 
in more early star formation in the modern semi-analytic 
models. At the same time, the use of the observed K-band 
z = luminosity function to normalize the PLE models 
has reduced the uncertainty due to dust and k-corrections 
in those models. The net effect is that the two scenarios di- 
verge significantly only at a higher redshift (z > 1.5) than 
predicted by KC98. Several other recent studies (Cimatti 
et al. 2002a; Firth et al. 2002; Kashikawa et al. 2003) 
have reached a similar conclusion. 

However, the morphologies of the observed objects in 
this redshift interval are inconsistent with the passive evo- 
lution hypothesis — the majority seem to be highly dis- 
turbed morphologically, and many are clearly interacting 
or merging (see also Daddi et al. 2003). Qualitatively, this 
is clearly more consistent with the hierarchical scenario. 
However, the quantitative disagreement between the num- 
ber of predicted and observed objects indicates that some 
ingredients in the models need to be modified, or else that 
some physics is missing. Further study of the nature of 
this population at 1.5 < z < 2.2, which forms a 'bridge' 
between the better-studied populations of 'normal' galax- 
ies at z < 1 and Lyman-break galaxies at z > 2.2, will 
certainly provide important new insights into some of the 
remaining mysteries of galaxy formation. 
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Fig. 2. — Observed frame color (R — K)ab vs redshift for the GOODS CDFS [left] and the semi-analytic mock catalog [right]. 
Tracks are shown for single age, stellar metallicity populations with ages of 13.5, 5.8, 3.2, 1, 0.5, and 0.1 Gyr (unreddened), from 
top to bottom. Large star symbols show the colors of spectroscopically confirmed objects from Daddi et al. (2003). The color 
distribution of galaxies with K < 22 in the redshift range 1.7 < z < 2.5 is shown in the right panel as histograms (diagonal 
hatched: GOODS; horizontal fill: mock catalog). 
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